Approximate Killing Fields as an Eigenvalue Problem 



Christopher Beetle 
Department of Physics, Florida Atlantic University, Boca Raton, Florida, 33431 

(Dated: July 13, 2007) 

Approximate Killing vector fields are expected to help define physically meaningful spins for non- 
symmetric black holes in general relativity. However, it is not obvious how such fields should be 
defined geometrically. This paper relates a definition suggested recently by Cook and Whiting to 
an older proposal by Matzner, which seems to have been overlooked in the recent literature. It also 
describes how to calculate approximate Killing fields based on these proposals using an efficient 
scheme that could be of immediate practical use in numerical relativity. 



oo : 
o . 

o ■ 
<^ : 



o 

I 



> 

00 

o 

oo 

o 



X 



Spacetime symmetries are essential for defining physi- 
cally important conserved quantities such as energy and 
angular momentum in general relativity. For example, 
when a vacuum spacetime admits a rotational symmetry 
generated by a Killing vector field , then a Komar-type 
integral P, over an arbitrary 2-sphere in spacetime can 
reproduce the physically well-defined angular momentum 
measured at infinity. When spacetime is axi-symmetric, 
but not vacuum, the difference between these integrals 
for a pair of different 2-spheres is precisely the ordinary 
angular momentum computed from the stress-energy of 
matter in the intervening space. If spacetime is not axi- 
symmetric, however, such formulae become rather am- 
biguous. They depend not only on the 2-sphere S over 
which one integrates, but also on the vector field used 
to define the integrand. 

These difficulties can be partially avoided under phys- 
ically favorable conditions. For instance, angular mo- 
mentum is well-defined at infinity in asymptotically flat 
spacetimes (see Q for a recent review), or on an appro- 
priate isolated horizon [1, modeling an isolated, non- 
dynamical black hole in a spacetime that may describe 
interesting dynamics in other regions. Essentially, these 
techniques identify preferred 2-spheres S (infinity, hori- 
zon, etc.) over which to integrate, and thereby reduce 
the ambiguity in defining the angular momentum. The 
resulting quasi-local formulae have the general Brown- 
York [6] form 



(1) 



where 5 is a (perhaps preferred) 2-sphere, Kat is the 
extrinsic curvature of a spatial slice E containing it, c/j" 
is a vector field on S, and dS'' is the area element of S 
within S. The basic problem remains: the vector field 
is arbitrary unless 5* has an intrinsic symmetry that 
can be used to select it. (Now, at least, that symmetry 
need not extend into the bulk of spacetime.) 

The horizons of black holes resulting from numerical 
simulations of astrophysical processes generally have no 
symmetry of any kind and therefore, seemingly, no pre- 
ferred vector field ip"". The problem is not that such 
surfaces have no reasonable definition for the angular 
momentum, but rather that they have infinitely many. 



There is one for every vector field tp"" tangent to the hori- 
zon. What is needed is a technique to pick a preferred 
vector field, and the obvious thing to do is to seek a ip"- 
that, in some sense, is as close as possible to a Killing 
field, even if none is present. This leads intuitively to the 
idea of an approximate Killing field. 

Motivated by the general issues discussed above, sev- 
eral groups have recently proposed elegant definitions of 
approximate Killing fields on 2-spheres. These include 
schemes based on Killing transport ^7] , conformal Killing 
vectors J3] , and most recently a minimization scheme by 
Cook and Whiting 9]. This paper revives an older ap- 
proach [13] due to Matzner based on solving an eigen- 
value problem, ft also suggests a novel adaptation of 
Matzner's approach to the specific problem of computing 
a preferred angular momentum for black holes in numer- 
ical relativity, and elucidates the intimate relationship 
between this scheme and that of Cook and Whiting. 

Let us begin with a brief review of Matzner's defini- 
tion of an approximate Killing field on a compact 
manifold M of dimension n equipped with a Riemannian 
metric gab- A continuous symmetry of gab is generated 
by a Killing vector field satisfying 



j!f^gab = 2V(^a^b) =0 



(2) 



throughout Af , where denotes the Lie derivative along 
and Va is the unique torsion-free connection on AI 
determined by gab- Taking a divergence, we see that any 
geometry with a continuous symmetry admits at least 
one non-trivial solution to the equation 



2VfcV('T^ = 0. 



(3) 



In principle, even when the geometry is not symmetric, 
we are still free to seek solutions for this equation. But 
generically we will not find any. 
Consider the eigenvalue problem 



(4) 



on a generic geometry. The operator Ak appearing here 
arises naturally in the transverse decomposition of sym- 
metric tensor fields on Riemannian manifolds [l^. A re- 
lated operator, denoted A^, arises in the same way from 
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the conformal KiUing equation, and plays a similar role 
in the transverse- iroce/ess decomposition of such fields. 
Its application to the initial-data pro blem in general rel- 
ativity is very well known indeed [l3| . 

Eq. ([4]) of course admits solutions only for a cer- 
tain spectrum of eigenvalues k, and zero may or may not 
be among these. Matzner establishes the following 
four results: the spectrum of eigenvalues k of Eq. ^ 
on a compact manifold is (a) discrete, (b) non-negative, 
(c) corresponds to a complete set of real vector eigen- 
fields and (d) contains k = if and only if the cor- 
responding eigenfield is a genuine Killing field. That 
is, the zero eigenspace of Eq. Q is precisely the finite- 
dimensional vector space of Killing fields of the metric gat 
on M. Therefore, Eq. ^ admits no solution if the met- 
ric gab on M has no continuous symmetries, as claimed 
above. However, we assert that the best approximation 
to a Killing field on a manifold with no actual symmetry 
is the unique vector eigenfield of Eq. with the min- 
imum eigenvalue k > 0. This is Matzner's definition of 
an approximate Killing field, and it has several desirable 
features. It exists generically, reduces to the correct an- 
swer when symmetries do exist, and, like a true Killing 
field on a symmetric manifold, is naturally defined only 
up to an overall (i.e., constant over M) scaling. 

Like any eigenvalue problem, Eq. (j4|) admits a varia- 
tional formulation. Recall the natural inner product 

(CO := / ClCe (5) 
Jm 

on the space of (complex) vector fields over M. Here, e 
denotes the canonical n-form volume element induced on 
M by the metric gab- We minimize the quadratic form 

QkK;^) := i (C - ^'^{{LO - 1), (6) 

where k is constant over M and here plays the role of 
a Lagrange multiplier. Minimizing this functional pro- 
duces the Euler-Lagrange equations 

A^r = «r and (C0 = 1, (7) 

the solutions of which are clearly the vector eigenfields of 
Eq. ([5]), normalized to unity in Hilbert space. 

Many variational problems are solved by initially solv- 
ing the first, differential equation in Eq. ([7]) for as 
a function of an arbitrary Lagrange multiplier k, and 
then using that result in the second, algebraic equation 
to impose the constraint and determine n. This does not 
happen for Eq. ([7]) because the differential equation is 
linear in and therefore leaves the overall scaling of 

undetermined. The second equation serves only to fix 
that scaling, and cannot also determine k. The Lagrange 
multiplier therefore must be fixed when we solve the first 
equation; for general k, no solution exists. This is hardly 
surprising since of course only true eigenvalues k allow 



us to solve Eq. ([4]) for However, it does make an ap- 
proach to Matzner's eigenvalue problem via a variational 
principle like Eq. ^ rather complicated. There is no al- 
gebraic equation to determine the Lagrange multiplier. 
Indeed, k is determined in this problem precisely by the 
condition that it be an eigenvalue of Ak, and there is no 
algebraic equation giving these. Minimizing Qk[^', i^) in 
Eq. ^ by solving the associated Euler-Lagrange equa- 
tions is neither easier nor harder than solving the eigen- 
value problem in Eq. 

Cook and Whiting's recent definition Q of an approx- 
imate Killing field uses a variational principle based on 
a quadratic form closely related to that of Eq. How- 
ever, it differs in a two important details. First, it focuses 
on the case where M ~ S* is topologically a 2-sphere, and 
restricts to be area-preserving: 

i?eea6 = (Vcf)eah = 0. (8) 

The motivation for this restriction arises from the tech- 
nical details of an eventual application to calculating the 
angular momentum of a non-symmetric black hole 
Second, it is based on a non-standard inner product 

(c,Ofl:=(c,i?o = £c:r^e- (9) 

These choices change the form, but not the basic content, 
of the resulting equations. They still describe a sort of 
self- adjoint eigenvalue problem. 

To restrict to area-preserving vector fields, it is easiest 
simply to recall that any divergenceless vector field 
on a 2-sphere topology is described by a unique scalar 
potential O such that 

t ^ {*d&T -.^ -t^^VbQ and ^ee = 0. (10) 

Now consider the restricted eigenvalue problem 

Kk^ ■.^{PAKP)t = kt, (11) 

where P denotes the orthogonal projection onto the sub- 
space of area-preserving vector fields within the Hilbert 
space of Eq. (O. A given = (*d6)" solves this equa- 
tion if and only if, for all other — (>i=d<I>)°, we have 

(*d$,AA'*de) K(*d$,*de}. (12) 

Integrating both sides by parts, and using positivity of 
the standard inner product of scalar functions over S, 
we find that Eq. (fTT|) is equivalent to 

A/fO := A2e-HV''(i?Vae) = KAe, (13) 

where A :— — V Va denotes the standard scalar Lapla- 
cian. We have shown that the scalar functions Q solv- 
ing Eq. (fT5)) generate, via Eq. PH)) . solutions of the 
restricted eigenvalue problem of Eq. pTjl . Because the 
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projection P does not typically commute with Aif, these 
do not generally solve Eq. (|4]) , and the restricted eigen- 
values k are generally distinct from the eigenvalues k in 
the full Hilbert space. In fact, we should generally ex- 
pect that Kmin > Kmin- Howcvcr, the area-preserving vec- 
tor eigenfield corresponding to this minimum restricted 
eigenvalue can also be considered a best approximation 
to a Killing field, albeit within a restricted class. 

To recover the Cook-Whiting approximate Killing 
field, we must repeat the previous calculation in the inner 
product of Eq. ©. The operator H^k is then no longer 
self-adjoint, but R^^ times Aif is. Accordingly, we seek 
vector fields = (^dO)" satisfying 



Kd*, R^^ Ak *de}^ = KR (*d$, *de)j 



(14) 



for all C° — (^d^)"^. Integrating by parts, and once again 
invoking positivity of the standard inner product of scalar 
functions, we find that Eq. is equivalent to 



(15) 



Although the notation here differs slightly, this is pre- 
cisely the Euler-Lagrange equation that Cook and Whit- 
ing find 0] by minimizing a quadratic form similar to 
Eq. dS]). Once again, the solutions (f^, kr) of this eigen- 
value problem generally differ from the solutions (^",k) 
of Eq. (m and from the solutions {C, k) of Eq. (fTT|) . 

Let us now make two technical comments. First, any 
constant function Q = c will give zero on both sides of 
Eqs. (fTT|) and (fT5| for all values of k ot kr, respectively. 
These are spurious solutions. They arise only because 
we have used potentials to describe the subspace of area- 
preserving vector fields. These solutions are ruled out by 
the the second condition in Eq. (fTO]) . which makes the 
correspondence between and Q an isomorphism. 

Second, the Cook-Whiting inner product in Eq ^ 
looks a little odd, but it is not immediately clear whether 
there is anything technically wrong with it. There cer- 
tainly can be problems. Recall that the scalar curvature 
in two dimensions varies as 



S^R : 



1 2 
2 
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under a perturbation 6gab of the metric. If this pertur- 
bation varies sufficiently rapidly over S, then the first 
term here can easily dominate the second, as well as the 
unperturbed, background value ^R. The result is that a 
generic spherical geometry, even if perturbatively close 
to a round sphere in the sense that dgab has small ampli- 
tude, can have regions of negative scalar curvature. (This 
is intuitively obvious if we imagine "pinching" the surface 
of a round sphere to create a small, saddle-shaped region 
of negative curvature.) On such geometries, the "inner 
product" of Eq. ([9]) is not positive-definite, and does not 
define a Hilbert space. However, in the space of all spher- 
ical geometries, there should be some finite region of "suf- 
ficiently smooth" perturbations of the round sphere for 



which the total scalar curvature remains everywhere pos- 
itive. In this region, there is no obvious problem with the 
Cook- Whiting scheme, but nothing particular to recom- 
mend it either. The question could presumably be settled 
[m by comparing qualitative features of the approximate 
KiUing fields computed from Eqs. pT|) and 

Matzner's eigenvalue definition of an approximate 
Killing field is unambiguous, universally applicable, and 
reproduces the usual Killing fields on a symmetric man- 
ifold. But it is not necessarily efficient in practice. In- 
deed, it would prohibitively expensive to solve any of the 
eigenvalue problems in Eqs. ([4]), (fTTj) or (fTSjl on the ap- 
parent horizon at every moment of time of a black hole 
in a numerical simulation. Fortunately, however, there 
is a simple approximation to speed the process up on a 
generic geometry. This approximation is based on the 
Rayleigh-Ritz method [15| , and works so long as we only 
want to find the lowest eigenvalue and the corresponding 
vector eigenfield. 

Consider the Rayleigh-Ritz functional 

F[(\ ^^'^^^ = -^^i (17) 
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on the full Hilbert space of Eq. ([5]), with the zero vec- 
tor removed. The local extrema of Eq. (|17p occur when 
f is a vector eigenfield of Ak, and the value of F[£] 
at each such extremum is the corresponding eigenvalue 
K. Note that the numerator here, which arises via in- 
tegration by parts of the second-order operator Ak in 
Eq. H]), is precisely one half the square integral of gab 
from Eq. ([2]). Thus, among all vector fields with fixed 
L^-norm on S, diffeomorphisms along the approximate 
Killing field modify the metric least in a quantifiable, 
sense. 

It is still not practicable to find the genuine abso- 
lute minimum of Eq. (fT7|) on the computer, which of 
course would yield Matzner's approximate Killing field. 
But one can approximate that minimum by minimizing 
F[£] within an appropriate space of trial vector fields. 
This idea is familiar from elementary quantum mechan- 
ics, where just such a variational principle is routinely 
used to approximate the ground-state wave-function of a 
complicated system. Unless the subspace of trial fields 
one chooses happens to be orthogonal, or nearly so, to 
the true minimum ^^^q of in all of Hilbert space, the 
dominant component of the minimizing trial field C^riai 
should lie along ^"j.^^, in Hilbert space. Most randomly- 
chosen trial spaces will not be orthogonal to Crruc- This 
idea allows us to approximate Matzner's approximate 
KiUing field. 

There is a natural candidate for the trial space of vec- 
tor fields in which to minimize Eq. (|17|) in the specific 
case M ~ S* of a 2-sphere horizon of a quiescent black 
hole in numerical relativity. One striking feature of many 
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recent numerical simulations (e.g, lj|) is that the hori- 
zons at late times often look fairly regular in the fiducial 
spacetime coordinates used to do the evolution. There- 
fore, it is natural to try a space of trial fields i^" based 
simply on those coordinates. A specific proposal follows. 

Use the fiducial spacetime coordinates in which the nu- 
merical evolution occurs to induce spherical coordinates 
(9, (f) on the black-hole horizon in some more-or-less nat- 
ural, but fundamentally ad hoc, way. Then, take the 
space of scalar trial potentials 
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where Yim{d,(f>) are the ordinary scalar spherical har- 
monic functions on a round sphere, 0'™ are arbitrary 
constants, and /max is a chosen cut-off. Each of these 
potentials generates an area-preserving vector field via 
Eq. (fTO|) . and this will be our trial space within the 
full Hilbert space of Eq. ([5]). Therefore, minimize 



F[Q] 



{*de,AK*de) _ {e,AKQ) 
(*de,*de) ~ (e,Ae) 



(19) 
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within the trial space of potentials given by Eq. 
Generally, we should expect that the minimizing po- 
tential will generate a vector field ^"^.j^^j fairly close to 
the minimum-eigenvalue area-preserving vector eigen- 
field ^^'j-uo of Eq. (fTTj) . This, in turn, should approxi- 
mate Matzner's approximate Killing field from Eq. 
To check the approximation, one could imagine increas- 
ing Zmax until doesn't vary much with the cut-off. 
Equivalently, one could use a fairly large cut-off — perhaps 
/max = 5 would bc cnough — from the start, and check 
that the amplitudes 0'™ are small for large /. If one 
prefers to approximate the Cook- Whiting approximate 
Killing field, one need only insert a factor of the scalar 
curvature R between the gradients in the denominator of 
Eq. (HID. 

There is one significant issue that has not been ad- 
dressed in this discussion. Even once an approximate 
Killing field ^° has been found from the eigenvalue ap- 
proach, it is still determined only up to overall normal- 
ization on S. For a proper rotational Killing field on 
a symmetric apparent horizon, the correct normalization 
would demand that the affine length of each Killing orbit 
should be 27r. It is not immediately clear what conven- 
tion might be used in the general case, without symme- 
try, to fix a normalization that goes over to this correct 
one in the limit of a symmetric manifold. This issue will 
be discussed more thoroughly, in the context of practical 
applications, in a forthcoming paper 11|. 
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The space of potentials in Eq. (|18|l is usually not orthog- 
onal, in the sense of Eq. (|10|l . to the space of constant 
functions on S. However, the key point is that standard 
properties of the ordinary spherical harmonics show that 
this space of trial potentials contains no actual constant 
functions. This is why we have taken /min = 1 in Eq. (|18|l . 
Therefore, Eq. (|10p maps our space of trial potentials 
faithfully to a space of trial vector fields with the same 
dimension, /max (/max -|- 2). 



